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To apply the scattering approach for the problem of AC transport through 
coherent quantum conductors, various partial density of states must be evaluated. 
If the global partial density of states (GPDOS) is calculated externally using the 
energy derivatives of the scattering matrix, the results are not precise unless 
the conductor has a large scattering volume. We propose a local formula for 
GPDOS which is suitable for any finite scattering volume. We apply this formula 
to compute the emittance of a two-dimensional quantum wire under the multi- 
mode and finite temperature condition. 
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The dynamic conductance of a quantum coherent mesoscopic system under 
a time dependent external field is the subject of recent interests]!], |2], |3|, [|, ^ ||. 
In contrast to dc-transport in the linear regime, where the internal potential dis- 
tribution inside a sample does not appear explicitly, the AC-response depends 
sensitively on the internal potential distribution. This internal potential is due 
to the charge distribution generated by the applied external AC-field at the leads 
and it has to be determined self-consistently|I]]. A particular useful approach to 
investigate AC transport properties of coherent quantum conductors, is to study 
the AC-response of the system to an external perturbation which prescribes the 
potentials in the reservoirs only|0, |]]. The external potentials effectively deter- 
mine the chemical potential of the reservoirs and the potential distribution in 
the conductor must be considered a part of the response which is to be calcu- 
lated self-consistently. In this approach, Biittiker and his coworkers [|l], |8| have 
formulated a current conserving formalism for the low frequency admittance of 
mesoscopic conductors. 

In the theory of Biittiker, Pretre and Thomas|l|, it is necessary to consider the 
Coulomb interactions between the many charges inside the sample, in order to 
preserve the current conservation. For a multi-probe conductor the low frequency 
admittance is found to have the form|| || G a p(u) = G a p(0) — iuE a/3 + 0(uj 2 ), 
where G a p(ti) is the dc- conductance, E a p is the emittance||, and a (or (3) labels 
the probe. The emittance E a p describes the current response at probe a due 
to a variation of the electro-chemical potential at probe j3 to leading order with 
respect to frequency lu. It can be written as|J 

E a/3 = d £ - D a p, (1) 

where the term dN a p/dE is the global partial density of states (GPDOS)|T0| 
which is related to the scattering matrix. It describes the density of states of 
carriers injected in probe (3 reaching probe a and is due to the response to the 
external perturbation. The term D a p is due to the Coulomb interaction of elec- 
trons inside the sample and is the response to the internal potential. D a p can be 
computed from the local density of states [p], || which is related to the electron 
dwell times through the relation J2 a D a /3 = Tdp/h where Tap is the dwell time 
for particles coming from the probe f3. Electric current conservation, namely 
J2 a G a p{oj) = 0, means that J2 a Ea@ — or equivalent ly|l], |TTf 

dNp _ x - dN a/3 x ^ . . 
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where dNp/dE is the DOS for electron coming from the probe (3. Clearly the 
current conservation is established since one realizes that ^ a dN a p/ dE is the 
physical quantity called injectance which is identical || to J2aD a p. 

The physical meaning and the important role played by the various partial 
density of states (PDOS) is the subjects of extensive discussions |TIJ. While in 
one-dimensional (ID) systems the PDOS can be evaluated analytically via the 
help of scattering Green's function, in 2D one is usually forced to use numerical 
methods due to the complexity of the problem 0, except for very special and 
exactly solvable cases [E^]. In the AC transport formalism outlined above, the 
GPDOS can be expressed approximately in terms of the energy derivative of the 
scattering matrix elements ||13||: 



dN, 



a/3 



dE 



— fdE(-f)Tr 
Am J v dE' 



ds, 



all 



^ dE 



dsj/3 
dE 



(3) 



In Eq.(^) f(E) is Fermi distribution function and s a p is the scattering matrix of 
dimensions M a x Mp where M a is the number of quantum channels supported by 
the lead a. Because for a given system one may be able to obtain the scattering 
matrix, Eq.(|3|) thus provides a practical means of computing the GPDOS. 

Using explicit 2D examples, it has been numerically demonstrated in Ref. 
[jjfl and analytically shown in Ref. ||, that the DOS dNp/dE computed from 
Eq.(|3]) is only accurate up to correction terms which are exponentially decaying 
functions of the scattering volume away from thresholds of successive transport 
channels. At the thresholds, the correction terms diverge. Hence precise current 
conservation can not be obtained from the explicit calculations [§, || |14| until the 
correction terms are added. A more serious problem is that we only know how to 
correct the DOS since we can compare the left and right hand sides of Eq. (Q), 
but we don't know how to distribute the correction terms among the PDOS of Eq. 
(||). The whole issue comes about since the external response in the theoretical 
formalism is really formulated with a time dependent perturbation at minus or 
plus infinity which are asymptotically far from the scattering region, while for 
every practical calculation one wants to consider what happens in some finite 



scattering volume pq . 

The purpose of the present work is to formulate a procedure which allows 
a precise determination of the GPDOS dN a p/dE for any finite scattering vol- 
ume. We have used a general approach similar to that employed by Christen 



and Buttikerjlql in computing the non-linear transport coefficients by eliminat- 
ing the GPDOS. This is possible if we assume that electric current conservation 
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is satisfied. In doing so all the required quantities for the emittance become local, 
thus can be calculated precisely using the new formula (see below) within a finite 
scattering volume. We have applied this formula to compute the emittance of a 
T-shaped quantum wire under the multi-mode and finite temperature condition. 
To proceed we recall from our previous investigations || H] of 2D quantum 



conductors, that the discrepancy between the DOS as computed externally or 
internally is related to the mode mixing of the 2D scattering. In 2D situations 
complicated mode mixing takes place. This mode mixing generates evanescent 
modes which can not propagate in the leads. For a scattering volume with a finite 
linear size L, the evanescent modes may "leak" out of the volume. However when 
we calculate the GPDOS from the scattering matrix using Eq. (|]), these "leaked" 
evanescent modes are not explicitly included, leading to a slightly inaccurate 
calculation. On the other hand, when we compute the local density of states 
internally (see below) using the scattering wavef unctions, all the modes, including 
the evanescent modes, are included. Indeed, as emphasized by Buttiker||, the AC 
transport formalism guarantees electric current conservation when the scattering 
volume is large enough to ensure that there is no electric field lines penetrating the 
surface of the volume. This condition is certainly violated due to the "leaked" 
evanescent modes when the volume is small. Since evanescent modes do not 
contribute to electric current (but does to the DOS), it seems to be natural to 
use the conservation law to eliminate the need of computing GPDOS externally. 
This will be our approach. 

Using the electric current conservation relation of Eq. (0) and the injectance 
formula! 



of J2 a Dap, we have 
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(4) 



This is consistent with the relationship between the GPDOS and local PDOS 
(LPDOS) 



dN^ 
dE 



dn a/ 3(r) 
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dn 
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(5) 



(6) 



is the LPDOS |fT0|- Hence the GPDOS can be computed locally through Eq.@ if 
we can obtain LPDOS. In general, for 2D systems it is very difficult if not impos- 
sible to obtain LPDOS using Eq. (||) since a complicated functional derivative 
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must be evaluated. In 2D this functional derivative can only be computed for 
specially simple systems ]|T2] . 



For ID systems, on the other hand, as shown in Ref. |HJ simplification to 
the LPDOS formula Eq. ([]) can be obtained (see Ref.|10| for details) via the 
Fisher-Lee relation^] between the scattering matrix and the retarded Green's 
function: 

Sq/3 = -8 aj3 + ih^/v a vpG(x a , xp) (7) 
where x a is the boundary of the scattering region. The functional derivative of 



the Green's function 5G/5U is given by [10 

8G(x a ,xp) 



5U(x) 



G(x a ,x)G(x,xp) 



Furthermore one can prove that for ID systems the following relation is true [10], 

G(x\,x)G(x,X2) — G(xi,X2)G(x,x) (9) 

for x\ < x < x 2 . Using Eqs. (0) - ®, it is not difficult to derive ]TIJ, for ID 
systems, the following expressions for the LPDOS 



dn. 



for a 7^ j3 and 



dn 



(313{ 
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(10) 



(11) 



where T Q/ g is the transmission coefficient from lead ft to a. In these results, the 
local DOS dnp{r)/dE is defined as 



drtt. 



dE 



E 



dn. 



af3{ 



dE 



(12) 



which is called the injectivity and it measures the additional local charge density 
brought into the sample at point r by the oscillating chemical potential at probe /3. 
In general, the injectivity can be expressed in terms of the scattering wavefunction 
aslH 



dn 



0{ 



dE 



hv 0n 



(13) 



where vp n is the velocity of carriers at the Fermi energy at mode n in probe (3. 
A related quantity, dn^^/dE, called emissivity, describes the local density of 
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states of carriers at point r which are emitted by the conductor at probe a. It is 
defined as 

dn*(T) _ x - dn aP {v) 
dE y dE ' 1 } 

It has been shown 0] that in the absence of a magnetic field the injectivity is equal 
to the emissivity. In the presence of a magnetic field the microreversibility of the 
scattering matrix implies that the emissivity into contact a in magnetic field B 
is equal to the injectivity of contact a if the magnetic field is reversed |J, 

dngjr, B) _ dn a (r, —B) 

dE ~ dE • (i5j 
Finally, dn(r)/dE is the local DOS given by 

dn(r) dn af3 (r) 

-dE- = ^^E- ■ (16) 



From these results, one is able to calculate the LPDOS via Eq. ([H]) and (|TT 
using the scattering wavef unctions. 

For 2D systems, the Fisher-Lee relation is of the form[p"8f 



(yp)dy a dyp (17) 

where Xan is the transverse wavefunction in lead a. In 2D, the equations similar 
to Eq. ([T0| ) and (|LTD do not seem to apply. This is because Eq. ([5j) does not hold in 



2D|19|]. However, since in a general 2D case one can not obtain analytical expres- 
sions for the scattering matrix and numerical computation are usually needed, it 
is thus enough to have a numerical prescription for calculating the LPDOS. Our 
numerical method makes use of a mathematical identity 



dV J 6U(r) ' v ' 

where ds a p/dV is calculated as follows[^(| ||, [12| : adding a constant potential 
V in the scattering volume and computing the scattering matrix formally to get 
s a/3 = S a/3(V), then taking the derivative and putting V = 0. Numerically the 
derivative can be easily carried out using finite differencing. 

Using Eqs.fll8|), (J5|), and (|J), the GPDOS is completely expressed by local 
quantities determined inside the scattering volume thus can be computed accu- 
rately for any system sizes, 



dN, 



a/3 



dE Alii 



1 fdE{-g)Tr 
7ii J v dE' 



t ds a/ 3 ^ ds^p 

Sal3 ~dV dv SaP 



(19) 



6 



This equation is also valid in ID. In the following, we shall compute the emittance 
from Eq. ([I]) using the Eq. (|I9"D for GPDOS. To obtain the quantity D a/3 we shall 
use the Thomas- Fermi approximation, in which case D a ^ is easily calculable [1], |8|] 
from Eqs. (0), (0), (0), and §L§, 



f 3 (dr^{T)/dE){dnp{T)ldE) 
D ^~J d r dn(r)/dE • (20) 

To illustrate the numerical procedure, we have computed the emittance of a 
T-junction quantum wire under the multi-mode and finite temperature condi- 
tions. The same system has been examined before for the single mode and zero 
temperature situation^ using externally computed GPDOS (i.e. Eq. (|3|)), thus 
the results here also provide an useful comparison. As shown in Fig. (1), the 
wire has two probes extending to x = ±oo while the scattering region is provided 
by the T-junction as shown being bounded by the two dotted lines. We assume 
that the boundaries of this ballistic conductor are hard walls, i.e. the potential 
V = oo. Inside the conductor the potential is zero. From now on we set h — 1 
and m = 1/2 to fix our units. 

Fig. (2) shows the emittance En versus incoming energy E in the first sub- 
band at zero temperature, where the GPDOS is computed using Eq.flT^) (solid 
line) or using Eq.(||) (dotted line). The difference in the two curves comes solely 
from the difference in GPDOS. Notice that near the second subband the diver- 
gence in the dotted curve is removed by using the locally computed GPDOS 
Eq. ([l9D as shown by the solid curve. As discussed previously [|J and shown in 
Fig. (2), there is also a divergence near E — E\ where E\ = n 2 is the first 
subband energy, if using the externally computed GPDOS. This problem is also 
overcome using Eq. (19) as shown by the solid curve. In fact since all quantities 
are computed locally if using Eq. (|19|), there is no need for any correction terms 
to the DOS. In Fig. (3a) we plot the emittance En up to the 2nd subband energy 
for three different temperatures T = 0.005i?i (dotted line), 0.05-Ei (dashed line), 
and O.lEi (solid line). The inset shows the same quantity up to 3rd subband. We 
have also depicted the dwell time r<ji (dashed line) in Fig. (3b) together with the 
DC conductance (solid line) and the emittance En at zero temperature (dotted 
line). There are eight resonant states located near the peaks of the dwell time. 
The behavior of the emittance En near the resonant energies £5 (the 5th one) and 
Sq (the 6th) which are basically inductive are different from that near the other 
resonant energies which are capacitive. This behavior is useful in explaining the 
resonance behavior in the transmission coefficient near resonant energies £5 and 
Sq. For instance, from the dwell time we know that near £ 5 there is a resonant 
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state. However, it is not clear whether it is a resonant transmission or resonant 
reflection since the DC conductance has a sharp peak and a dip very close to 
each other. The fact that the emittance near £ 5 is inductive indicates that it is 
a resonant transmission. Finally, as the temperature is turned on, the peaks or 
valleys in the emittance curve diminish gradually. When the temperature reaches 
T = 0.2Ei the interference pattern is washed out completely. This is expected 
since a finite temperature tends to smear out the quantum resonances. 

In summary, we have proposed a numerical procedure and formula for com- 
puting the global partial density of states which is precise for any finite scattering 
volume of a quantum conductor. As GPDOS plays a most important role in the 
AC transport theory, our result provides a useful tool for further numerical in- 
vestigations of the dynamic admittance. In this formulation the electric current 
conservation is satisfied automatically. This formulation, especially Eq . (|l8|) , also 



applies to the non-linear transport jig, PI - Applying the new procedure to a 
T-shaped junction, the divergences of the emittance at each subband edge are 
removed. 
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Figure Captions 

Figure 1. Schematic plot of the quantum wire system: The wire width, the side-sub 
width, and height are fixed at W = 1. The two dotted lines separate the 
scattering region from the two probes. 

Figure 2. The emittance En versus incoming energy E in the first subband at zero 
temperature, where the GPDOS is calculated using Eq.(|T9"D (solid line) or 
using Eq.(||) (dotted line). 

Figure 3. Figure (3a): the emittance En versus energy for three difference tempera- 
tures T = 0.005-Ei (dotted line), 0.05.&L (dashed line), and O.lEx (solid line), 
up to the second subband. Inset, the same quantity up to the third sub- 
band. Figure (3b): the conductance Gn (solid line), the dwell time (dashed 
line), and the emittance En at zero temperature (dotted line) where the 
dwell time and the emittance have been scaled by a factor of eight to fit in 
one figure. 
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